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I.  INTRODUCTION 


The  particular  problem  of  estimating  gun  barrel  temperature  profiles  in 
the  vicinity  of  high  energy  gas  flow  has  been  examined  experimentally  by  local 
investigators,  and  the  general  problem  of  mathematically  constrained,  234 
unsteady,  high  temperature  heat  flux  is  treated  in  the  open  literature*  ’  ’ 
Historically,  the  preliminary  layout  of  liquid  propellant  rocket  combustion 
chambers  has  often  employed  the  Schmidt  diagram3’  to  estimate  the  critical 
metal  temperatures  in  the  structure.  This  is  developed  as  a  finite  difference 
approximation  of  the  second  order,  first  degree  partial  differential  equation 
of  physics  in  heat  flow  context  and  is  usually  presented  graphically.  The 
graphical  analysis  is  an  alternative  to  the  LaPlace  Transform  treatment'  which 
is  unwieldy  in  execution.  It  also  competes  with  the  differential  analyzer 
technique.  Rigorous  solutions  to  the  second  order  axisymmetric  problem  in  the 
form  of  mapped  finite  element  presentations  of  the  detailed  transient  flow 
phenomena  are  promised**  for  the  proximate  future. 

A  specific  and  comprehensive  solution  to  the  barrel  problem  has  recently 
been  delivered  in  the  form  of  mathematical  presentation  by  Polk,^  and  this  in 
the  explicitly  usable  form  of  a  moderate  level  computer  program  (Hewlett- 
Packard  9845)  for  direct  application.  The  Polk  treatise  is  not  currently  in 
BRL  report  format,  however,  and  the  program  is  still  being  refined.  A 

T.  L.  Brosseau,  "An  Experimental  Method  for  Accurately  Determining 
Temperature  Distribution  and  the  Heat  Transferred  in  Gun  Barrels »  BRL-R- 
1740,  September  1974.  AD  HB000171L. 

2  Mark  W.  Zemansky,  Heat  and  Thermodynamics,  McGraw-Hill  Book  Company  Inc., 

New  York,  1957. 

3  Max  Jacob,  Heat  Transfer,  Vol.  1 ,  John  Wiley  and  Sons,  New  York ,  1949. 

4  E.  F.  Quigley,  "One  Dimensional  Transient  Temperature  and  Stress 
Distribution  Produced  in  0.375-  and  0.500-  inch  Thick  7075A1-T6  Flat 
Plates  by  Fourteen  Nuclear  Thermal  Environment ,"  BRL-MR-2173 ,  April 
1972.  AD  *901995. 

0  Frank  Kreith,  Principles  of  Heat  Transfer,  International  Text  Book 
Company,  Scranton,  PA.,  1963. 

ft 

George  P.  Sutton,  Rocket  Propulsion  Elements,  John  Wiley  and  Sons,  New 
York,  1956. 

Theodore  V.  Harman  and  Maurice  A.  Biot,  Mathematical  Methods  in 
Engineering ,  McGraw-Hill  Book  Company  Inc.,  New  York  1940. 

9  P •  L.  Versteegen  and  F.  D.  Varcolik ,  "Heat  Transfer  Studies  in  Gun  Tubes," 

ARBRl-CR-00393  (Science  Applications  Inc.,  McLean,  VA .)  March  1979.  AD  M069649 

9  J.  F.  Polk,  "An  Algorithm  for  Heat  Transfer  in  Gun  Barrels,"  Transactions 
of  the  Twenty-Fifth  Conference  of  Army  Mathematicians,  ARC  Report  80-1, 

1980 . 
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complete  FORTRAN  program,  based  on  a  Calspan  reference,^  is  also  available 
but  has  yet  to  be  introduced  into  the  local  Control  Data  Corporation  facility. 

In  direct  graphical  exposition,  the  Schmidt  plot  consists  of  a  large 
scale  drafting  wherein  the  results  are  recovered  as  measured  values  rather 
than  discrete  numbers.  This  report  offers  the  details  of  a  desk  top 
calculator  (Hewlett-Packard  97)  procedure  algebraically  transposed  from  the 
finite  difference  scheme  to  calculate  heat  transfer  coefficients  in  a 
ballistics  environment  from  an  overdetermined  set  of  equations.  The  data  is 
from  experiment, and  an  implicit  assumption  is  that  the  driving  gas  temperature 
is  constant  over  the  selected  time  interval  of  the  calculation.  The 
appropriate  driving  gas  temperature  can  then  be  determined  by  regressive 
substitution. 

Thermodynamic  units  always  deserve  a  special  note.  The  examples  are 
presented  in  the  International  System  of  Units,  the  British, and  in  normalized 
parameters.  The  British  units  are  employed  simply  because  the  materials 
properties  values  are  most  widely  available  in  this  system, and  it  therefore 
offers  direct  and  convenient  application.  The  SI  is  proposed  as  a  standard, 
and  the  normalized  borrows  from  the  Jacob-Hawkins-3  practice  with  an  additional 
modification  to  eliminate  geometric  scale  factors. 

II.  PROCEDURE 


A.  General  Statement 

The  differential  equation  describing  one  dimensional,  unsteady, 
conductive  heat  flow  through  a  solid  is 


I  il 

a  3t 


3  T 

— j  (rectangular  slab). 

3x 


(1) 


where  a  is  thermal  diffusivity, 
T  is  temperature, 
t  is  time, 
x  is  distance. 


Approximate  solutions  to  the  differential  equations  can  be  obtained  by  solving 
the  finite  difference  analogue. 

,  AT  A  2T 
It  x 

a  At  2  (2) 

Ax 


10  "CYCOND"  Program ,  and  "HTC"  Program ,  oommunioation  via  Mr.  Kovacs ,  DRDAR- 
SEy  Picatinny  Arsenal. 


10 


with 


At  equal  to  the  finite  difference  in  time, 

2 

Ax  equal  to  the  finite  difference  in  distance, 

A^T  equal  to  the  time  variable  effecting  a  change  in  T,  and 
A^T  equal  to  the  distance  variable  effecting  a  change  in  T. 

By  then  defining  "t"  as  the  number  of  time  increments  (At),  and  "n”  as  the 
distance  increment,  the  direct  application  form  of  the  equation  becomes 


.  T  (t+1)  -  T  (t)  T,  Aj.(t)  -  2  T  (t)  +  T,  ,.(t) 
1  n _ n  _  (n+1 ) _ n _ (n-1 ) 


a  At  ,2 

Ax 


(3) 


Note:  Throughout  this  report,  the  notation  T  t  and  TR(t)  will  be  used 
interchangeably.  Parenthesis  will  be  used  only  when  added  clarity  is 
necessary. 

Appendix  A  transcribes  the  finite  differences  to  algebraic  formulation  and 
illustrates  the  Schmidt  plot  by  numerical  example. 


In  cylindrical  coordinates  the  finite  difference  equation  becomes 


.  T  (t+1)  -  T  (t) 
1  __n _ n 

a  At 


(T(n+I)(t)  *  T(n-l)(t>  ~  2  Tn(t)} 


Aj 


(4) 


where  r  =  e^  , 

a  A  •  _  Al¬ 
and  Aj  -  — —  , 


for  any  radius  "r"  with  "j"  defined  in  terras  of  "r". 

Appendix  B  presents  the  cylindrical  transcription  and  a  numerical  example  for 
illustration. 


B.  Application  to  a  Gun  Barrel 

The  gun  barrel  is  represented  as  a  long  thick-walled  annular  cylinder, 
Figure  1,  with  a  suddenly  imposed  high  temperature  gas  flowing  through  the 
bore.  To  estimate  the  wall  temperature  gradient,  three  assumptions  are 
required : 


A  gas  temperature  history. 

Univalued  properties  of  the  medium. 

Radiation  and  convection  accountable  by  the  modification  of  the  value 
for  unit  surface  conductance. 


where 

r^  is  the  inner  cylinder  radius 

r  is  the  outer  cylinder  radius 
o 


Figure  1.  Schematic  of  Infinitely  Long  Cylinder 


C.  Wall  Surface  Temperature 


The  procedure  used  to  develop  the  equations  which  specify  the  temperature 
of  the  face  of  a  slab  after  exposure  to  a  high  temperature  gas  consists  of: 

( 

1.  Construction  of  the  Schmidt  plot,  which  is  geometrically  defined  in 
terras  of  selected  gas  and  metal  properties;  i.e.,  the  heat  transfer 
coefficient  (h),  the  reference  temperature  (T^),  the  conductivities 
of  the  gas  and  the  metal  (k),and  the  diffusivity  of  the  metal  (a). 

2.  Transformation  from  the  geometric  to  the  equivalent  algebraic 
expression. 

In  this  context,  as  shown  in  Figure  2. 


AT  ■  T  o  -  T  o. 
o  o  n 


m 


Ax 

Ax+Ax  * 
o 


or 


(5) 


_ 1 _ 

1  +  (Ax  /Ax) 

o 


(6) 


where  Ax 

o 


k 

—  by  construction  of  the  Schmidt  plot  and 
h 


Ax 


(2aAt) 


(7) 


where  At 


duration  of  heat  exposure 
number  of  time  intervals 


For  the  gun  problem,  the  duration  of  the  heat  exposure  may  be  determined  from 
the  in-bore  projectile  trajectory  history  and  an  adiabatic  flame  temperature 
established  from  a  definition  of  the  propellant  composition.  These  are 
tenuous  criteria,  of  course,  but  are  within  the  existing  "state  of  the  art." 
From  Figure  2  and  Table  1,  the  temperatures  at  the  wall  (Plane  1)  are 
algebraically  fixed  for  the  first  ten  time  intervals  in  terras  of  equivalent 
Schmidt  plot  distances  corresponding  to  the  specific  thermal  resistances  on 
each  side  of  the  wall  plane.  Table  1,  obtained  from  Figures  3  through  7, 
presents  the  algebraic  resolution  of  the  temperature  profile  for  the 
interstitial  planes.  Table  2  gives  the  corresponding  temperatures  for  the 
second  plane  -  based  on  the  same  reasoning.  These  particular  formulae  neglect 
the  influence  of  curvature  since  the  wall  penetration  (Ax)  approaches  zero. 


Figure  2.  Wall  Surface  Temperature  Gradient 


Figure  4.  Graphical  Determination  of  T  3 


mAT0 


e  6.  Graphical  Determination  of  T.. 7 


U  Determination  of  T,9 


1  *1  '-^1  IS.1  W}  ^  IK  "K  ^TT^fTT^rr^TT  ■-".  «.  ■■'  V  yff  ■;".’  V  7  ■■■  V  <■  «■■  ».l 


TABLE  1.  ALGEBRAIC  SPECIFICATION  OF  TjD 


Ti1 


Tl3 


Tl5 


Ti7 


Tj9  - 


where 


Ti°  *  V 

Tj2  -  m  A  T0  +  TjO 

Tj4  *  V2  ®  A  T0(3-m)  +  TjO 

Tj6  *  V4  m  ATo(7.5-4.5m  ♦  m2)  +  T^o 

Tj8  ■  mATo(17.5-14.5m  ♦  6m3  -  m3)  +  T^o 

T, 10  -  vr  mAT  (39. 375-40. 625m  +  23.75m2-7.5m3+m4)+T1o 
1  10  o  1 

<t>  is  as  defined  above, 
n 


■  m  AT  4>,+T1o 
01  1 

-  m  AT  4»-+T.o 

o  3  1 

•  m  AT  $_+T,o 

o  5  1 

*  m  AT  ♦,+T,o 

o  7  1 

-  m  AT  4>q+T,o 

o  9  1 


TABLE  2.  ALGEBRAIC  SPECIFICATION  OP  T2n 


T22  -  T23  *  V2  m  ATq  +  TjO 

T24  *  T25  "  *4  m  AT0(3»5_,nJ  +  T2o 
T26  -  T27  -  j  mATQ(9.5-5m  +  m2)  +  T2o 

T2®  *  T2^  *  T6  m^^0(24.25  -  19.75  m  +  6.5  m2-m3)  +  T 

where  6  is  as  defined  above, 
n 


■  m  AT  0-+T-O 

o  3  2 

*  1  AT  Qr+T-o 

0  3  Z 

■  m  AT  8  +T«o 

o  /  z 

*  »  4T„VV 
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The  calculating  technique  is  demonstrated  by  example,  using  the  HP-97,  in 
Appendix  C. 

Evaluation  of  h  from  Experimental  Data 

If  an  experimental  wall  temperature  distribution^*4  can  be  established, 
it  is  possible  to  read  the  temperature  rise  at  any  two  times  and  solve  the 
algebraic  set  of  equations  of  Table  1  to  eliminate  Tq.  The  resulting 
expressions  are  single  valued  in  "ra"  and  solvable.  Using  the  most  convenient 
divisor , 


V  -  V 

raAT  <}>, 
o  1 


Vo  (3-m)  =  Cm)  , 


(8) 


v  -v 

mAT 

o  1 


=  l/2  (7.5-4.5m+ro2)  = 


4>5(m)  , 


(9) 


Ti7  -  V 
mAVi 


*  ^  (17.5-14.5m+6m2-m3)  =  <^(ra)  , 


(10) 


Tl9  -  Tl0 
mAT 

o  1 


jy  ( 39. 375-40. 625m+23 . 75m2-7 .5  m3  +  m4)  =  <J>9(m),(ll) 


Ttn  -  TjO 

whereby  — ^ -  =  m  ^(ra)^.  (12) 

o 

The  solution  iterates  "ra"  until  the  experimental  and  calculated  T^n  differ  by 
less  than  a  pre-selected  error  bound.  The  heat  transfer  coefficient,  h  ,  is 
then  deduced  from  the  definition  of  "ra." 


Appendix  D  contains  a  sample  calculation  and  printout  of  the  HP-97  program  for 
illustration. 

Determination  of  Tgn  from  Experimental  Data 

Plainly,  the  extension  of  the  procedure  is  to  examine  the  average  value 
of  the  driving  gas  temperature  over  the  subintervals,  i.e.,  those  within  the 
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elapsed  time  Interval.  Using  values  of  T^l  and  Tj3,  the  Intersection  of  the 
rays  by  the  construction  of  Figure  4  locates  the  value  of  TQn  for  the 
subinterval  between  time  1  and  time  3.  The  remaining  intersections  are 
similarly  determined  except  for  the  period  between  time  0  and  time  1.  T^l  Is 
directly  determined  by  the  physical  constraints  of  the  finite  difference1 
equation.  Figures  8  through  12  illustrate  and  Tables  3  and  4,  based  on  the 
following  identities,  summarize  the  procedure.  Table  3  considers  m  -  constant 
while  Table  4  allows  a  variable  m. 


v 

—  -  m,  0<t<l  ,  (13) 


Tj3  -  T23 

V  "  T23 


m  ,  l<t<3 


(14) 


TjS  -  T25 

T  n  -  T5 
o  2 


m  ,  3<t<5  , 


Tj7  -  T27 

V  '  V 


-  m  ,  5<t<7  , 


Tj9  -  T29 

T  n  -  T9 
o  2 


m  ,  7<t<9  . 


(15) 


(16) 


(17) 


T2n  is  listed  in  Table  2  and  Tn0  commonly  adds  to  zero. 


For  the  subinterval  temperature  rise  and  a  corresponding  "m,"  the  local  heat 
transfer  coefficient  can  be  calculated  by  the  methods  demonstrated  in  Appendix 


TABLE  3 


Time  Interval 


From  To 


TABLE  4.  ALGEBRAIC  SPECIFICATION  OF  T0n  FOR  ■  +  CONSTANT 


III.  RESULTS 


The  data  from  References  11,  12  is  used  to  examine  some  of  the 
implications  of  the  equations.  Figure  13x^  provides  the  temperature  history 
of  the  interior  barrel  surface  of  a  bench  mounted  20  mm  weapon  firing  TB-1 
propellant.  The  thermocouple  was  located  immediately  forward  of  the 
chamber.  For  the  idealized  conditions  previously  specified,  i.e.,  constant  TQ 
equal  to  the  adiabatic  flame  temperature  of  the  propellant  and  univalued  h  and 
with  the  properties  given  in  Table  D-l,  the  heat  transfer  coefficient 
corresponding  to  a  maximum  wall  surface  temperature  rise  of  842 *F  (450 *C)  is 
found  to  be  10,953  B/hr  ft^F  (1656  cal/hr  cm^  C).  This  result  is  obtained  by 
iterating  "m"  and  T^9  -  T^o  until  the  calculated  agrees  with  the  experimental 
temperatures.  Figure  14  compares  three  different  driving  gas  temperatures  for 
the  cases  of  an  assumed  constant  heat  transfer  coefficient  and  adiabatic  flame 
temperature  source  and  the  alternate  cases  where  the  heat  transfer 
coefficient  is  allowed  to  vary  within  the  subintervals — in  one  instance  being 
an  averaged  constant  over  the  full  interval  and  in  the  other  being  coupled  to 
the  local  effective  gas  temperature.  Figure  15  shows  the  corresponding  heat 
transfer  coefficients  for  this  example. 

IV.  CONCLUSIONS 

With  reliable  digital  temperature  recording  to  one  millisecond  resolution 
and  current  desk  top  calculating  facility,  the  Schmidt  plot  approach  can  be 
used  to  establish  the  effective  heat  transfer  coefficient  for  experimental 
data  where  the  transport  properties  of  the  system  may  be  assumed  constant. 


Artillery  Ammunition  Master  Calibration  Chart ,  Material  Testing 
Directorate  Report  137 S,  15th  Revision ,  Aberdeen  Proving  Ground  t  MD 
1973. 
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Unpublished  Test  Datat  Interior  Ballistics  Division t  Ballistic  Research 
Laboratory  APG,  MD  1980. 


Experimental  Wall  Temperature  History  for  Special  Propellant  Round 


TIME  INTERVAL  SEQUENCE 


Figure  14.  Driving  Gas  Temperature  Profiles 
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APPENDIX  A 


TEMPERATURE  GRADIENT  IN  A  TWO  DIMENSIONAL  PLANE  WALL 

For  the  simple  geometry  represented  by  a  semi-infinite  flat  plate 
exposed  to  a  suddenly  imposed  constant  temperature  gas  on  one  surface, 
the  Schmidt  plot  takes  the  form  given  on  Figure  A-l.  The  initial  distance 
increment,  Ax_,  is  made  proportional  to  the  film  coefficient  conductance 
while  the  following  nine.  Ax,  characterize  the  metal.  Within  the  given 
calculation,  the  number  of  distance  increments  is  determined  by  the 
storage  capacity  of  the  calculator.  There  is  no  limit  to  the  number  of 
time  cycles. 


T2°jj\  '3,-3 


INTERFACE  PLANE  LOCATIONS 


Figure  A-l.  Schmidt  Plot  for  Schematic  Slab  Temperature  Profile 


The  problem  is  to  establish  the  interface  time-temperatures  in 
order  to  insure  that  the  wall  material  strength  criteria  is  not  exceeded. 
Table  A-l,  following,  illustrates  the  stepping  progression  in  the  finite 
difference  parameters  and  Figure  A-2  presents  the  results  of  a  specific 
example. 


TABLE  A-l.  STEPPING  PROGRESSION  IN  TEMPERATURE  TIME  PARAMETERS  FOR  SLAB 
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There  ;ire  two  branches  to  the  algorithm.  One  branch  describes  the 
conditions  at  odd-numbered  time  cycles  such  that 


Tjt  =  i^t-l  + 


4X  (T0-T2t-l) 


Ax  +  Ax 


(A-l ) 


T2t  =  T2t-1 


(A-2) 


T  t 
3 


T-t-1  +  T.t-1 

2  4 

2 


(A-3) 


T  t  =  T.t-1 

4  4 


v  - 


T.t-1  +  T  t - 1 

4  6 


T  t  =  T.t-1 

6  o 


T.t-1  ♦  Tut-1 

T  t  =  — _ - _ 

7 1  2 


V  -  V-1 


(A-4) 

(A-5) 

CA-6) 

(A- 7) 

(A-8) 


T  t 
9 


r8t_1  +  1 10t_1 
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(A-9) 


T 


10 


V-1 


(A-10) 
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The  second  branch  describes  conditions  at  even-numbered  time  cycles  such 
that 


V 

-  v-1 

(A— 11) 

'  t 

Tjt-1 

+  V-1 

2 

(A- 12) 

3* 

=  v-1 

(A- 13) 

4t 

V-1 

+  V-1 

(A-14) 

2 

S* 

=  Tst-1 

(A-15) 

1  f 

V-1 

+  T?t-1 

(A-16) 

b 

2 

7Z 

-  I7t-1 

(A- 17) 

T,t-1  +  T,,t  - 1 


The  magnitude  of  the  time  cycles  is  determined  from  the  initial  definitions 
such  that 


Axi_ 

2aAt 


1 


(A-21) 


or 


At  » 


» 


where 


a  = 


k_ 

pc 


(A-22) 


and  k  is  thermal  conductivity, 

c  is  specific  heat,  and 
p  is  material  density. 


Whereby 


At 


cAx 

2pk 


(A-23) 


Example 


The  calculating  technique  is  demonstrated  by  example.  Assume  a  semi-infinite 
steel  wall  12  inches  (30.2  cm)  thick  suddenly  exposed  to  a  5000°F  (2760°C) 
gas  flow  over  one  face.  After  a  measured  interval  of  26  sec,  the  gas 
temperature  is  reduced  to  ambient.  Table  A-2  defines  the  initial  conditions 
and  operating  parameters  for  the  program  input.  Figure  A-2  presents  a  Schmidt 
plot  of  the  results.  The  HP-97  listing  follow  Figure  A-2. 


TABU  A- OLHNIIION  OF  BOUSUARV  CONDITIONS  FOR  SAMPLE  SLAB  PROBLEM 


sec  I  At  =  Ax^/2 


•  V*  •>  "V  r.’.'V. 1  v 


111-97  INITIAL  KKCISTKK  CONTENTS 

Registers  R  through  R()  enter  the  initial  temperatures  at  the  interstitial 
planes. 

R^  indicates  the  number  of  time  cycles  required  to  elapse  2.6  minutes 
(.for  this  example). 

R„  gives  the  distance  increment  fpr  the  steel. 

R,  gives  the  sum  of  the  steel  and  equivalent  film  increments. 

R„  and  R|.  indicate  the  initial  temperature  at  the  terminal  plane. 

U|  enters  the  number  of  time  increments  completed  (usually  0). 

R„  gives  the  effective  ambient  temperature  of  the  gas  during  the  cooling. 


OUTPUT 


The  printed  output  presents  the  temperature  in  the  enumerated  sequence 
of  the  interstitial  planes.  The  terminal  plane  temperature  is  given  in 
the  I)  register  and  the  number  of  time  cycles  completed  in  the  1  register. 
The  real  time  elapsed  is  this  number  times  the  time  interval. 


SAMPLE  INPUT  AND  OUTPUT 


The  program  entry  employed  with  the  HP-97  is  given  below, 
consistant  system  of  units  is  suitable. 


Primary 

Registers 

R 


Initial  Values 

To 

T1 

r. 


T, 


,VI 
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R, 
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‘8 

1 9 
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Ax 
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10 
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AMBIENT 


Example 

2730  N 
280  N 
280  N 
280  N 
280  N 
280  N 
280  N 
280  N 
280  N 
280  N 
6 

.208  CAL 
1.208  CAL 
208  N 
280  N 
0 


280  N 


Output  after  15  time  cycles 
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T. 


280  N 
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—  30.19  cm - 


Figure  A-2.  Schmidt  Plot  for  Sample  Slab  Probl 


PROGRAM  LISTING 
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APPENDIX  B 


TEMPERATURE  GRADIENT  IN  A  LONG  THICK  WALLED  CYLINDER 

In  cylindrical  coordinates  the  suddenly  imposed  temperature  is 
shown  in  Schmidt  plot  presentation  as  Figure  B-l,  which  is  a  transfor¬ 
mation  from  the  physical  plane  to  that  of  the  operating  parameters.  Fro 
the  equations  of  I I -A,  the  transformed  equivalent  film  distance  is  equal 

to  — =-  for  equal  radial  increments  transposed  to  logarithmically  spaced 
rh 

increments  in  Aj . 


oc 
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iu 
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Figure  B-l.  Schmidt  Plot  Schematic  for  Cylindrical  Wall  Temperature 
Profile 

The  problem  is  that  of  the  slab  to  establish  the  time-temperature 
distance  gradient.  Table  B-l  illustrates  the  stepping  procedure  in 
finite  difference  parameters,  and  Figure  B-l  plots  the  results  of  a 
specific  example. 


E  B-l.  STEPPING  PROGRESSION  IN  TEMPERATURE  TIME  PARAMETERS  FOR  CYLINDER 


The  odd-numbered  time  cycles  are  described  by 


A  fixed  time  increment  is  found  such  that 


The  calculating  technique  is  demonstrated  by  example.  Assume  an  infinitely 
long  stainless  steel  cylinder  1.0  inches  (2.54  cm)  inside  diameter  and 
2.5  inches  (6.35  cm)  outside  diameter  suddenly  exposed  to  a  constant 
high  temperature  gas  flow  through  the  bore  for  1.31  sec.  The  properties 
of  the  materials  and  the  critical  transport  properties  are  given  in 
Table  B-2  with  the  HP-97  listing  immediately  following. 


TABLE  B-2 . 
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HP-97  REGISTER  CONTENTS 

The  geometric  radial  boundaries  are  specified  and  equal  radial 
increments  are  automatically  converted  to  logarithmic  parameters  within 
the  program.  The  corresponding  initial  temperatures  must  be  entered  as 
well  as  an  elapsed  time  span. 

Registers  Rg0  through  Rg^  enter  the  initial  temperatures  at  the  radial 
interstitial  planes. 

RA  indicates  the  number  of  time  cycles  required  to  elapse  seconds  (for 
this  example) . 

Rg  gives  the  inner  radius. 

R  gives  the  outer  radius. 

R^  indicates  the  effective  ambient  temperature  of  the  gas  during  the 
cooling  phase. 

Rg  enters  the  initial  temperature  of  the  terminal  plane. 

Rj.  is  given  the  number  1  for  the  initial  calculation.  The  program 

automatically  converts  this  figure  to  that  of  the  number  of  time  cycles 
completed  at  any  stage  in  the  calculation. 

Note  that  the  program  starts  with  the  secondary  register  filled  and  the 
primary  register  open. 


OUTPUT 

The  printed  output  presents  the  temperature  in  the  enumerated 
sequence  of  the  radial  interstitial  planes.  The  terminal  plane 
temperature  is  given  in  the  D  register  and  the  number  of  time  cycles 
completed  in  the  I  register.  The  real  time  elapsed  is  this  number  times 
the  time  interval. 
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PLANE  WALL  SURFACE  TEMPERATURE 


APPENDIX  C 


PLANE  WALL  SURFACE  TEMPERATURE 

Under  the  idealized  conditions  as  previously  specified,  the  wall 
temperature  progression  is  established  by  algebraic  determination  as 
expressed  in  Table  1.  For  example,  assuming  an  adiabatic  flame  tempera¬ 
ture  of  M-2  propellant  as  the  driving  gas  medium,  and  estimating  the 
time  required  for  the  projectile  to  clear  the  barrel  as  .003  sec,  and 
a  heat  transfer  coefficient  of  1845  B/hr  ft^  F,  with  the  physicals  as 
given  in  Table  C-l,  the  wall  temperature  at  the  five  time  intervals  can 
be  read  directly  from  the  printout.  The  complete  HP-97  listing  follows 
on  p. 


INITIAL  REGISTER  CONTENTS 

This  program  requires  a  constant  driving  gas  temperature  T  ,  an  initial 
wall  temperature,  the  thermal  conductance  and  diffusivity  of  the  wall 
Material,  and  a  convective  heat  transfer  coefficient  corresponding  to  an 
assumed  elapsed  time  increment.  It  presumes  a  slab  wall  and  calculates 
the  wall  temperature  growth  from  an  initial  slab  depth  of  approximately 
.040  in  (1  mm) . 

R.  indicates  the  driving  gas  temperature,  T  n. 

Rg  gives  the  initial  ambient  temperature,  Tno. 

R^  is  the  elapsed  time  in  hours. 

R2  enters  10,  the  number  of  time  increments  for  this  example. 

R3  is  the  average  convective  heat  transfer  coefficient  over  this 
range . 

R^  gives  the  metal  thermal  conductivity. 

R^  is  the  metal  thermal  diffusivity. 

OUTPUT 


The  printed  output  gives  the  wall  temperature  for  ten  sequential 
increments  of  time  over  the  elapsed  time  interval. 


TABLE  C-l.  DEFINITION  OF  BOUNDARY  CONDITIONS  FOR  SAMPLE  PLANE  WALL  PROBLEM 


The  program  entry  employed  with  the  HP-97  is  given  below.  Any  consistant 
system  of  units  is  suitable. 


Primary 

Initial 

Registers 

Values 

Example 

5596  F 

R 

T  n 

a 

0 

rb 

T  o 

70  F 

R1 

n 

t 

-7 

8.33x10  hr 

R2 

e 

10 

10 

R3 

h 

1845  B/hr  ft 

R4 

k 

16  B/hr  ft 

R5 

a 

.49  ft2/hr 

Output 

after 

3-millisecond  time 

V 

= 

v 

246  F 

Tl3 

= 

V 

331  F 

Tl5 

= 

Tl6 

394  F 

V 

= 

Tx8 

446  F 

^9 

= 

T^O 

490  F 

67 
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PROGRAM  LISTING 

891 

X 

-35 

041 

* 

-55 

892 

8 

82 

842 

PRTX 

-14 

893 

i 

n  4 
it 

843 

RCL2 

36  02 

094 

RCLB 

36  13 

844 

*2 

53 

895 

♦ 

-55 

045 

STGS 

35  09 

896 

PRTX 

-14 

046 

RCL2 

j6  82 

897 

RCL9 

36  03 

047 

4 

04 

898 

X* 

53 

001 

*LBL8 

21  li 

848 

• 

-62 

899 

8 

00 

002 

PRE6 

16-13 

049 

5 

05 

188 

X 

-35 

083 

RCL1 

36  01 

858 

X 

-35 

181 

ST07 

w5  0, 

804 

RCL2 

36  82 

851 

CHS 

-33 

182 

RCL9 

36  65 

885 

4 

'tat 

852 

RCL  3 

36  03 

103 

CHS 

fc*. 

886 

37  CC 

35  13 

853 

♦ 

-55 

184 

6 

06 

887 

RCL5 

36  05 

854 

? 

07 

18  5 

6 

00 

888 

X 

-35 

055 

• 

• 

-62 

186 

X 

-35 

889 

■t 

£ 

02 

056 

5 

05 

187 

RCL  7 

Jo  Ot 

818 

X 

-35 

05  7 

* 

-55 

188 

+ 

_  err 
wo 

811 

IX 

51 

858 

RCL4 

36  04 

189 

S  TO? 

?r  r.  ? 
WW  if  1 

812 

STGl 

35  01 

05  9 

X 

-35 

110 

RCL9 

36  05 

813 

RCL4 

36  84 

068 

4 

04 

111 

1 

01 

814 

RCL  3 

36  83 

061 

■i 

-34 

112 

9 

03 

815 

• 

-24 

862 

RCLB 

36  12 

113 

8 

00 

816 

STUD 

35  14 

■  863 

+ 

.rr 

.. 

114 

X 

RCl7 

-35 

817 

RCL1 

36  81 

,  064 

PRTX 

.  . 

*T 

115 

36  6  7 

818 

-24 

065 

RCL9 

36  32 

116 

♦ 

-  rr 

wo 

819 

1 

01 

866 

RCL2 

36  02 

117 

S707 

35  67 

828 

+ 

-55 

867 

X 

-35 

118 

RCL2 

36  03 

821 

J/X 

CO 

ft/ft. 

068 

sm 

35  82 

119 

3 

1 73 

822 

STQ2 

35  32 

069 

CHS 

128 

2 

63 

823 

RCLR 

36  11 

070 

ST07 

35  07 

121 

5 

65 

824 

RClB 

36  12 

871 

RCL9 

36  05 

122 

X 

->>5 

825 

-45 

872 

6 

Sc 

123 

CHS 

.  0*0 
ta-ft» 

826 

S703 

35  03 

873 

X 

-35 

124 

RCL7 

Jo  3  i 

827 

RCL2 

36  02 

874 

* 

.rr 

wo 

125 

♦ 

-  C-r 

ta/ta* 

8  28 

X 

-35 

075 

sro  ? 

**5  0 1 

126 

$7  07 

35  C7 

829 

ST  04 

35  04 

876 

RCl2 

36  03 

127 

3 

03 

838 

RCLB 

36  12 

877 

1 

31 

128 

1 

31 

831 

♦ 

-55 

878 

4 

84 

129 

5 

35 

032 

PRTX 

-14 

879 

• 

-62 

138 

* 

_  rr 
ft/ 

833 

3 

03 

888 

5 

85 

131 

RCL  4 

3t  S-f 

834 

RU2 

36  32 

881 

X 

_  -r 

WO 

132 

X 

-35 

835 

-45 

882 

Cfi6 

n 

133 

A 

til 

836 

RCl4 

36  04 

883 

RCL7 

36  07 

134 

£ 

e: 

837 

X 

-7* 

wJ 

884 

+ 

.rr 

ftfft/ 

135 

i 

3£ 

838 

2 

00 

Ota 

885 

1 

■T  < 

t* 

136 

T 

"^t 

839 

•“4 

tat 

886 

7 

ti 

137 

RClB 

3o 

848 

RCLB 

*7  -•  4 

JG  ik 

887 

* 

“C»ta. 

138 

♦ 

_  rr 

%/ w 

088 

C 

J 

r 
v  . 

139 

PR  TX 

4  t 

089 

♦ 

_  rr 

148 

RTH 

ft.  t 

090 

RCL4 

f , 

r.  - 


68 
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APPENDIX  D 


HEAT  TRANSFER  COEFFICIENT  FROM  EXPERIMENTAL  DATA 

Given  the  conditions  recorded  in  Table  D-l,  the  HP-97  input/output 
and  the  complete  program  listing  follows.  This  program  contains  a 
recursive  feature  such  that  the  initial  estimate  of  h  is  used  to 
calculate  T^9,  which  is  then  compared  with  the  experimental  value, 
whereupon  a  recalculation  using  a  modified  value  of  h  is  used  to  reform 
T^9  until  it  is  within  the  acceptable  range  selected.  The  cyclic  converging 
values  of  T^9  are  printed  in  order,  and  the  result  presents  the  final  h  with 
the  corresponding  temperature  progression. 


DEFINITION  OF  BOUNDARY  CONDITIONS  FOR  EXPERIMENTAL  PROPELLANT  ROUND 


HP-97  INITIAL  REGISTER  CONTENTS 
Register  0  enters  the  range  of  temperature  convergence. 

Register  1  is  the  elapsed  time  in  hours. 

Register  2  is  10 

Register  3  allows  the  estimate  of  h,  the  film  coefficient. 

Register  4  gives  k,  the  conductivity  of  the  metal  wall. 

Register  5  indicates  a,  the  diffusivity  of  the  metal  wall. 

Registers  6-9  are  zero 

Register  A  is  the  driving  gas  temperature  T  n,  usually  taken  as  the 
local  adiabatic  flame  temperature. 

Register  B  is  the  initial  wall  temperature  T^o. 

Registers  6,  D  and  E  are  zero. 

Register  I  gives  the  wall  temperature  rise  in  “C.  If  the  calculation  is 
not  carried  out  in  British  units,  steps  153  to  159  in  the  listing  must 
be  adjusted  to  delete  the  conversion  of  T^n  from  ®F  to  °C. 

OUTPUT 

The  Intermediate  printout  is  the  temperature  rise  for  each 
permutation  of  h.  For  a  slowly  converging  operation,  the  program  can 
be  stopped  at  any  point  and  a  revised  h  entered  in  Register  3. 

The  final  printout  gives  the  complete  results  of  the  terminal 
calculation. 

Register  0  shows  h. 

Register  1  gives  T^l. 

Register  2  gives  T^3. 

Register  3  gives  T^5. 

Register  4  gives  T^7. 


Register  5  gives  1,9. 


Register  6  repeats  Register  I. 

Registers  7-9  are  zero. 

Register  A  repeats  the  driving  gas  temperature  TQn. 

Register  B  repeats  the  initial  wall  temperature  T^o. 

Registers  C  and  D  have  no  significance  to  the  result. 

Registers  E  and  I  present  the  final  calculated  and  the  reference 
temperature  rise  respectively. 
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The  actual  program  entry  is  given  below.  With  the  exception  noted  for 
Register  I,  any  consistant  system  of  units  is  suitable. 

Primary 


Regi sters 

Initial  Values 

Example 

R 

T 

2  C 

o 

R1 

z 

t 

e 

1.94xl0_7hr 

R_ 

Constant 

10 

2 

R3 

h 

12000  B/hr  ft2F 

R4 

R5 

k 

a 

16  B/hr  ft  F 
.49  ft2/hr 

R6 

0 

0  . 

R7 

0 

0 

*8 

0 

0 

R_ 

0 

0 

9 

ra 

V 

4398  F 

rb 

V 

70  F 

Rc 

0 

0 

rd 

0 

0 

re 

0 

0 

RI 

V  -  TjO 

450  C 

Intermediate 

values  of 

Tx9 

1  -  TjO  (°C) 

Output  after 

iteration. 

Ro 

h 

10953.2  B/hr  ft2F 

R! 

V 

410  F 

R2 

T  3 

567  F 

R3 

V 

678  F 

R4 

V 

767  F 

R5 

Tx9 

842  F 

R6 

Tj9  -  T  1  (calc) 

449.9  C 

R, 

0 

75 
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APPENDIX  E 


ALGEBRAIC  SPECIFICATION  OF  TQn 

The  program  labeled  E-3  refers  to  the  determination  of  the  driving  gas 
temperature  according  to  the  requirements  of  Table  3 .  Program  E-4  similarly 
refers  to  Table  4.  The  example  and  data  of  Appendix  D  are  used  for 
illustration. 

INITIAL  REGISTER  CONTENTS 

For  both  programs: 

R^  enters  the  wall  temperature  at  time  1 . 

R2  gives  the  wall  temperature  at  time  3 . 

R3  Indicates  the  wall  temperature  at  time  5. 

R4  is  the  wall  temperature  at  time  7 . 

R5  enters  the  wall  temperature  at  time  9. 

OUTPUT 

The  corresponding  driving  gas  temperatures  are  given  in  stated  sequence. 


f 


Primary 

Registers 

Initial 

Values 

Example 

R1 

T11 

410  F 

R2 

Tl3 

567  F 

R3 

Tl5 

678  F 

R4 

Tl7 

767  F 

R5 

Tx9 

842  F 

Output 

• 

Program  E-3 

To1 

820  F 

T03 

1751  F 

To5 

2056  F 

V 

2275  F 

V 

2430  F 

Program  E-4 

V 

820  F 

V 

1751  F 

To5 

2468  F 

V 

3278  F 

3469  F 


1  *.»  V 

*  v W  %  l  1  ^  v\  ^ 

as 

c& 

n 

S'/--. 

fry 

kV-", 

i 

881 

tLBLE 

to  A  AW 

882 

RCLl 

36  81 

’*,*7^  * 

883 

a 

to 

82 

884 

885 

X 

PRTX 

-35 

4  4 

“At 

H 

886 

RCL1 

36  61 

88? 

3 

83 

888 

X 

.?r 

ww 

889 

ST07 

?P  37 

WW  iM 

• .  ■* 

818 

RCL2 

Jo  &2 

011 

2 

02 

812 

X 

“wW 

.9, 

813 

RCL7 

JO  <)( 

. 

814 

- 

815 

CHS 

■22 

8 16 

ST08 

35  SO 

81? 

RCL1 

jo  01 

|g 

818 

819 

X* 

RCL8 

53 
36  60 

028 

‘24 

821 

PR  TX 

-14 

* . 

822 

RC-L1 

36  81 

823 

4 

84 

p 

824 

825 

X 

ST07 

-35 
35  87 

026 

RCL1 

36  61 

027 

RCL2 

36  02 

■» 

828 

9 

83 

.' v”  \ 

829 

X 

-35 

1  -*• 

838 

CHS 

“22 

831 

RCL7 

36  07 

832 

t 

“55 

833 

ST07 

35  07 

034 

RCL2 

36  02 

035 

X* 

r? 

WW 

036 

2 

02 

• 

83? 

X 

-  7r 

Ww 

838 

CHS 

- 

*  toto 

V  v- 

039 

RCL7 

36  07 

848 

♦ 

_  rr 
ww 

% 
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PROGRAM  LISTING  E-4 

089 

CHS 

‘toto 

341 

S70? 

35  07 

098 

RCL7 

JO  0l 

■  842 

RCL2 

jo  02 

091 

+ 

rr 

Ww 

843 

RCL3 

36  03 

892 

ST07 

Jw  0? 

844 

X 

-35 

893 

RCL2 

36  02 

045 

8 

00 

894 

RCL4 

36  04 

04  6 

X 

-35 

895 

X 

-35 

847 

RCL7 

36  07 

896 

S 

88 

848 

♦ 

_  r  r 
ww 

897 

X 

-35 

849 

ST00 

35  00 

098 

RCL 7 

JO  Ci 

850 

RCL1 

36  01 

899 

+ 

_  rr 

ww 

051 

3 

03 

188 

ST07 

Jw  Cl 

052 

X 

-35 

181 

RCL3 

36  S3 

853 

$707 

35  07 

182 

X'‘ 

C7 

ww 

854 

RCL2- 

36  82 

183 

8 

06 

855 

1 

01 

104 

X 

*»r 

ww 

056 

8 

00 

185 

CHS 

‘toJ 

85  7 

X 

-35 

186 

RCL7 

36  07 

858 

CHS 

-22 

187 

♦ 

_  e*r 

v  v 

859 

RCL7 

36  07 

108 

ST06 

35  OC 

060 

♦ 

-55 

189 

RCL4 

jo  04 

861 

ST07 

WW  dl 

118 

1 

01 

862 

RCL3 

jo  Cj 

111 

6 

06 

863 

8 

06 

112 

X 

_  7c 
ww 

864 

X 

-35 

113 

ST07 

Jw  0  > 

865 

RCL7 

36  07 

114 

RCL3 

36  03 

866 

* 

_  rr 
wW 

115 

1 

n  4 

WA 

867 

RCL8 

36  0O 

116 

6 

66 

868 

i 

-24 

117 

A 

-  7C* 
ww 

869 

UX 

52 

118 

CHS 

‘22 

878 

PRTX 

.  4  4 

it 

119 

RCL7 

36  07 

871 

RCLl 

36  01 

128 

♦ 

.rr 

Ww 

872 

RCL3 

36  03 

121 

RCLl 

36  61 

873 

X 

-35 

122 

+ 

-55 

874 

CHS 

to  A. 

123 

6707 

35  07 

875 

ST07 

35  0i ■ 

124 

RCL2 

JO  Cj 

876 

RCLl 

36  01 

125 

6 

06 

877 

RCL4 

36  04 

126 

A 

RCL? 

-35 

878 

X 

“ww 

127 

JO  iJ. 

879 

2 

82 

128 

+ 

.rr 

Ww 

888 

X 

-35 

129 

ST  07 

35  07 

881 

RCL7 

jo  tir 

138 

RisLj 

36  03 

882 

♦ 

ww 

131 

S 

02 

883 

ST07 

35  07 

132 

A 

__  7c 

V  w: 

884 

RCL2 

36  82 

133 

CHS 

toto 

885 

RCL3 

7 -t? 

Jv  Ww 

134 

RCL7 

7_-  :T 

wO  O. 

886 

X 

"  Jw 

135 

+ 

-rr 

ww 

887 

2 

02 

136 

RCL0 

36  00 

888 

X 

-  7r 
ww 

13? 

r 

■V  .4 

‘wt 

84 


139 

PRTX 

148 

RCL1 

141 

5 

142 

X 

143 

ST  07 

147 

RCL7 

148 

t 

149 

ST07 

150 

RCL3 

151 

1 

152 

6 

153 

X 

154 

RCL7 

155 

+ 

156 

ST07 

157 

RCL4 

163  RCL4 


166  ST07 

167  RC LI 


170  ST 08 

171  RCL2 

172  1 

173  6 

17  4  x 

175  RCL8 

176  + 

177  ST08 

178  RCL3 

179  6 

180  4 

181  x 

182  RCL8 


35  OS 

36  02 
01 


'Jw 

35  OS 

36  03 
06 
04 

-35 
36  OS 


184 

RCL  5 

36  05 

185 

X 

-35 

186 

RCL7 

36  0? 

187 

+ 

-55 

188 

STO0 

35  00 

189 

RCL  5 

36  35 

198 

6 

06 

191 

4 

04 

192 

X 

-35 

193 

ST07 

35  07 

194 

RCL2 

36  02 

195 

4 

04 

196 

X 

-35 

197 

RCL? 

ob  0i 

198 

+ 

~55 

199 

ST07 

35  0? 

280 

RCL3 

36  S3 

201 

2 

02 

282 

4 

04 

283 

X 

_  7S- 

wv 

284 

RCL? 

36  or 

205 

♦ 

■55 

286 

ST07 

35  07 

287 

RCL4 

36  04 

212 

RCL  7 

213 

♦ 

214 

2 

215 

X 

216 

RCL1 

217 

3 

218 

X 

219 

* 

228 

RCL8 

221 

T 

222 

i/X 

223 

PRTX 

224 

RTH 

APPENDIX  F 


UNITS  DESIGNATION 

For  the  given  examples,  the  properties  are  stated  in  the  British, 
S.I.  and  normalized  expression.  The  non-dimensional  geometric  "caliber" 
is  most  useful,  with  many  applications  in  fluid  dynamic  analysis,  and 
the  extension  to  a  normalized  mass  (weight)  has  been  a  convenience  in 
rationalizing  the  performance  of  grossly  different  flight  projectiles.^”! 
Such  manipulative  procedures  have  physical  limitations  which  are  defined 
in  the  framework  of  their  employment. 

The  normalization  parameters  are  given  below.  Dividing  weight 
(force)  quantities  by  the  specific  weight  of  water  converts  them  to 
length  quantities  and  the  mechanical  equivalent  of  heat  transforms  heat 
units  into  mechanical  units.  In  addition,  this  report  introduces  a 
common  temperature  ("Normalized)  based  on  the  absolute  scale. 


Length 


linear  dimension 


reference  dimension 


CAL 


Force 


force 


.  3 


specific  weight  of  water  x  ref  dim 


CAL' 


Mass 


mass 


.  3 


specific  mass  of  water  x  ref  dim 


CAL' 


Heat 


heat  x  mechanical  equivalent 
specific  weight  of  water  x  ref  dim^ 


CAL 


Temperature  “Normalized  =  .9  “Kelvin  =  .5“Rankine 


^”*W.F.  Donovan, 
Penetrators »" 


"One  Factor  Affecting  the  Dispersion  of  Long  Rod 
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LIST  OF  SYMBOLS 


a  thermal  diffusivity  of  metal  wall 
c  specific  heat  of  metal  wall 

e  exponential  constant 

e  elapsed  quantity  (when  appearing  as  a  subscript) 

h  film  heat  transfer  coefficient 

j  natural  logarithm  of  "r" 

k  thermal  conductivity  of  metal  wall 

Ax 

m  Ax  +  Ax 

0 

n  sequential  location  of  interstitial  plane  (appearing  as  subscript) 

r  radius 

t  time 

x  distance 

cal  calorie 

cm  centimeter 

ft  foot 

hr  hour 

sec  second 

B  British  Thermal  Unit 

C  Centigrade 

F  Fahrenheit 

L  Length 

N  Normalized 

T  Temperature 

CAL  Caliber 

A  Difference  indicator 

At  Finite  difference  in  time 

AT  Finite  difference  in  temperature 

Ax  Finite  difference  in  distance 

<j>  Operator  defined  as  used  in  text 

0  Operator  defined  as  used  in  text 


LIST  OF  SYMBOLS  (continued) 
p  Density  of  medium 

q  Heat  flux 

r.  Inner  radius 

1 

r  Outer  radius 

o 

t  Elapsed  time 

•f  H 

t  Time  at  n  interstitial  plane 

n 

t  Range  of  convergence  for  temperature  calculation 

T  t  Temperature  at  n  interstitial  plane  at  time  "t" 

Ax'  Initial  distance  increment  in  Schmidt  plot 

-o 

AT  Initial  driving  temperature  difference 

HP-97  Hewlett  Packard  Model  97  Calculator 

FORTRAN  Computer  language  acronym  for  formula  translation 

CDC  CONTROL  DATA  CORPORATION 

CALSPAN  Formerly  Cornell  Aeronautical  Laboratory 

SI  STANDARD  INTERNATIONAL 

BRL  BALLISTIC  RESEARCH  LABORATORY 
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